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We consider the one-dimensional target search process that involves an immobile target located 
at the origin and searchers performing independent Brownian motions starting at the initial 
positions x = (xi, X2, ■ ■ ■ , xn) all on the positive half space. The process stops when the target 
is first found by one of the searchers. We compute the probability distribution of the maximum 
distance m visited by the searchers till the stopping time and show that it has a power law tail: 
' Pjv(m|a;) ~ Bn{x\X2 ■ ■ .XN)/m^'^^ for large m. Thus all moments of m up to the order (TV — 1) 

I are finite, while the higher moments diverge. The prefactor Bm increases with faster than 

exponentially. Our solution gives the exit probability of a set of A'^ particles from a box [0, L\ 
through the left boundary. Incidentally, it also provides an exact solution of the Laplace's equation 
in an Ai'-dimensional hypercube with some prescribed boundary conditions. The analytical results 
are in excellent agreement with Monte Carlo simulations. 
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r-| ■ I. INTRODUCTION 

O , 

The probability distribution of the maximum of a single one-dimensional Brownian motion (and its variants such as a 
Brownian bridge or an excursion) over a fixed interval of time [0, t], has a long history in the probability literature [H-Q]. 
The statistics of the maximum has diverse applications. One example is the Kolmogorov-Smirnov test in statistics that 
^ \ is used to compare, in a nonparametric way, two different probability distributions [1, Q . Similarly, the distribution of 
C/3 . the global maximum of a discrete-time random flights (including Levy flights) has also been studied in the probability 
literature (lol . [flj , with more recent applications in computer science jl2j | , physics ^IS*! and chemistry [l4| . 

In statistical physics, there has been a recent revival of interest in related problems in the context of the distribution 
of the maximal height, measured with respect to a reference point, of (1 + l)-dimensional fluctuating interfaces [Tsl - 
[iTj . In the stationary state of a finite sample of size L, such fluctuating interfaces are often described by a Brownian 
^ • bridge in space over an interval [0,L], albeit with certain global constraints [l^. The statistics of maximum has 
Q [ also been computed for continuous-time subdiffusive processes p^ . IToj and has been used to analyse single particle 
O . trajectories |20| . 

The distribution of the maximum for a single Brownian motion (or its variants such as bridge, excursion etc.) 
, has been extended to many Brownian motions, including certain strongly interacting random walkers, e.g. non- 
^ • intersecting, so-called vicious random walkers l2lM24ll. (The latter problem has an intriguing connection to the 
\ Gaussian ensembles of the random matrix theory [2l|,[2J,[25.) For independent walkers, the results on the distribution 
\l of the maximum have recently been used to compute the mean perimeter and the mean area of the convex hull of N 
independent planar Brownian motions [2^ [27| . 
. These results on Brownian motion and its variants represent rare exact analytical results for the extreme value 
■ statistics of correZated random variables, a subject of increasing current interest [28|. However, all these results about 
,, the distribution of the maximum, for a single or multiple walkers, have been derived in the case when one considers 
the walkers over a fixed interval of time [0,t]. An interesting variation of this problem, with several applications, arises 
when the interval [0,t] is not fixed, but itself varies from realization to realization, i.e., one observes the walker (or 
walkers) over a time interval [0, tg] where the stopping time ts of the process itself is a random variable. For example, 
ts may represent the first-passage time (through the origin) of a walker. 

To be more precise, consider first a single Brownian walker that starts at time i = at position xi > 0. The position 
cd ' xi{t) of the walker evolves via the continuous-time stochastic equation, dxi/dt — ?7i(t) where 771 (t) is a Gaussian white 
noise with mean (771 (t)) = and a correlator {r]i{t)rji{t')) = 6{t — t'). The process stops at the stopping time tg when 
the walker hits the origin for the first time (see Fig. [1]). 

Let m be the maximum displacement of the particle till the stopping time tg . The statistics of the random variable 
m is interesting and it represents an example of the so called first-passage Brownian functional [1^ . The problem is a 
toy model of 'random search', where the origin represents a fixed 'target' and the Brownian walker represents a random 
searcher. The search is called off when the searcher finds its target and m represents the maximum distance travelled 
by the searcher before it finds its target. For concreteness, we shall mostly used terminology related to random 
search, althou gh there are several applications of this problem. For example, in the context of trapping [sol - ls^ or 
predator-prey [33| models the origin may represent an immobile target (prey) and the Brownian walker may represent 
a diffusing chemical trap (predator) . The stopping time ts is then the reaction time or the survival time of the prey 
and m denotes the maximum distance the predator travels before finding its prey. In the context of the directed 
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FIG. 1: The trajectory (red line) of a single Brownian walker starting initially at xi till the stopping time ts at which it hits 
the origin for the first time. The maximum distance travelled by the particle till ts is denoted by m. 



Abelian sandpile model in (1 + l)-diniensions [34| . m represents the maximum lateral size of an avalanche [35( 1 . The 
random variable m also plays an important role in characterizing the so called staircase polygons [s^ . In the context 
of queueing theory, where the position of the walker represents the length of a queue, m represents the maximum 
length of a queue during the so called busy period [35,, ji6\ . 

The probability density function (pdf) Pi{m\xi) of m (for fixed xi) can be easily computed (ssj and it turns out 
to be a pure power law: 

Pi{in\xi) = — ^; ■m>xi. (1) 

While this pdf is evidently normalized to unity, the average (m) and higher integer moments are infinite! The 
cumulative distribution of the maximum is given by 



Qi{L\xi) = Prob [m < L\xi] — I 



X\ 

Pi{m\xi) dm =1 — . (2) 

1^ 



This distribution has a very simple interpretation: it just represents the exit probability of a Brownian particle j37| . 
starting at < xi < L, from a box [0, L] through its left boundary at 0. 

In this paper, we study a generalization of this search problem where there is still one fixed target at the origin, but 
there are searchers who perform independent Brownian motions on the x > axis, starting at the initial positions 
X = (xi, X2, . . . , xn)- The position Xi{t) of the i-th walker evolves with time t via the Brownian evolution 

(3) 

where rji{t) is a Gaussian white noise satisfying {rji{t)) = and = ^(^ — t')^i,j- Since the walkers are 

independent they can cross each other. The process stops at a stopping time ts when the origin is hit for the first 
time by any one of the walkers (e.g., the second walker (red) in Fig. [2]). Note that tg varies from one history of the 
process to another. 

In the context of chemical kinetics (3l| where the problem is generally referred to as the 'target annihilation' 
problem, various generalizations of this problem have been investigated including e.g. the situation where the target 
itself diffuses |38l - l43 |. In the following, we shall limit ourselves to the case of immobile target and focus on the statistics 
of the maximum distance m (from the target) travelled by any of the walkers till the stopping time ts when the target 
is found. Thus m denotes the distance of the farthest point on the x axis visited by any one of the walkers till ts. 
Clearly m is a random variable fluctuating from one realization of the process to another. Our object of interest is the 
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FIG. 2: The trajectories of = 5 independent Brownian walkers starting at initial positions xi, X2, xs, Xi and x^ till the 
stopping time is when one of the walkers (the second one (red) in this figure) hits the origin. The maximum displacement 
along the x direction till ts (undergone, e.g., by the third particle (green)) is denoted by m. The cumulative probability 
Qn{L\^x) — Prob [m < L\x\ also represents the exit probability of the first particle from a box [0, L\ through its left boundary. 

probability density P/v(m|a;) of this maximum distance m, given the number TV of walkers and their initial positions 
X. Thus m provides an estimate (worst-case) of the distance that needs to be covered by a team of N walkers to find 
a fixed target. 

As in the single searcher case, let Qn{L\x) — Prob[TO < L\x\ = PN{m\x) dm be the cumulative probability that 
the maximum m till ts is less than or equal to L. This cumulative distribution of the maximum can be interpreted 
as the solution of a different problem as in the = 1 case. Consider, for instance, a slightly different problem where 
again we have a set of N independent walkers, but now inside a box [0, L], starting at the initial positions x. Let us 
define the exit probability as the probability that the first particle that exits the box [0, L] does so through (and not 
through the upper boundary at L), see Fig. [5] As in the N = 1 case, this exit probability is precisely the cumulative 
distribution Qn{L\x) of the maximum m till the stopping time in the semi-infinite system, as it counts all those 
events where one of the trajectories hits the lower boundary before hitting the upper boundary at L while all the 
others stay inside the box [0, L] till this event of first-hitting the origin. 

We will see that for this seemingly simple one-dimensional model of independent walkers, the statistics of m has a 
rich and nontrivial dependence on the number N of walkers. This is partly due to the fact that the same stopping 
time ts for all the walkers effectively introduces a correlation between the trajectories of the walkers, even though 
each executes an independent Brownian motion. While for A^ = 1 the solution is simple, it becomes rather nonrivial 
even for A = 2! 

Let us first summarize our main results. We compute the pdf PN{m\x) exactly for all A^ > 1 by a path counting 
(or path integral) method. We show that, for arbitrary A^ > 1, the pdf of the maximum has an asymptotic power-law 
tail 



where the prefactor Bn has a nontrivial A^ dependence which we compute explicitly. For A = 1, we have Bi — 1 and 
the asymptotic result in is actually valid exactly for all m> xi. For N — 2, we will see that 



PAr(m|x) ~ Bn 



XiX2 ■ . - Xn 



(4) 



as m — >■ c» 




(5) 



In particular, for large A', the prefactor Bn grows faster than exponentially 



BnO:^N - \n{N) 



(6) 



TT 
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Our asymptotic result ([4]) indicates that for A'' walkers, integer moments of m up to order (A'' — 1) are finite, while 
higher integer moments are infinite. Evidently, as N increases, the distribution becomes narrower and narrower as 
expected but it does so in a nontrivial fashion. 

Recently, the cumulative distribution of the maximum m till the first-passage time ts, or cquivalently the exit 
probability Q{L\x) from the box [0,L] through the origin, was studied [43] for a generic self-affinc stochastic process 
x{t) starting at the initial position x. The process x{t) typically grows with time as x{t) ^ where H is the Hurst 
exponent. This power law growth of distance with time makes the process self-afhne. An example is the ordinary 
Brownian motion where H = 1/2. For such a generic self-affine process, it was argued (isj that the cumulative 
distribution of the maximum Q{L\x) = Prob(rn < L\x) ^ 1 — A (x/L)'^ in the limit x/L — where A is a constant. 
The exponent (p was found to be related to the persistence exponent 9 via the scaling relation (p = 0/H [i^. The 
persistence exponent 9 characterizes the late time power law decay of the survival probability, i.e., the probability 
that the process stays on the positive half-axis up to time t 4J]. Thus the pdf of the maximum decays for large m 
as, P{m\x) ~ m~'^~^ with 4> = 9/ H . The exact result (|4]) shows that if we think of the assembly of N independent 
Brownian motions as a single self-afhne stochastic process in the TV-dimensional space, then (f) — N . We will see later 
that the persistence exponent for this collective process \s 9 = N /2 and the Hurst exponent H = 1/2. Thus our exact 
result for this model supports the general scaling relation (j) — 9/H found in (43i] . 

The paper is organized as follows. In Section 2, we provide a simple heuristic argument in favour of our main result 
(U]). This argument is not sufficient to compute the prefactor _Bjv exactly for all N. However, we show that this 
heuristic argument becomes asymptotically exact for large N and one can extract the limiting behavior oi Bn for 
large N using an extreme value argument. In Section 3, we set up the general method for computing the cumulative 
distribution Qn{L\x) of the maximum m. This requires solving Laplace's equation in an A^-dimensional space with 
appropriate boundary conditions. We present explicit solutions for the cases = 1 and N = 2. In Section 4, we 
present an alternative path counting method that is more general, physically transparent and provides explicit results 
for all A'^ > 1. In Section 5, we present numerical results to verify our analytical predictions. Finally in Section 6, 
we conclude with a summary and a list of interesting open problems. Some of the details of the computations are 
relegated to the appendices. 



II. HEURISTIC ARGUMENT 



We begin with a simple heuristic argument in favour of Consider the semi- infinite geometry with the indepen- 
dent Brownian motions, evolving via starting at initial positions {a;i,X2, . . . ,XAr}. Let tg be the stopping time 
when one of the walkers hits the origin. The probabality distribution of tg can be computed exactly as follows. As 
an input to solving the A^-particle problem, let us first consider a single Brownian motion starting initially at xq > 
and let p{xo,t) be the survival probability, i.e., the probability that the walker does not hit the origin up to time t. 
This can be easily computed by various standard methods and is given by the well known formula [29|, |45[ 

p{xo,t) = eri ( ^] ; erf(y) = ^ r e'^" du. (7) 



2tJ Vtt 







Turning to A^ walkers we note that since the walkers are independent, the probability that none of them hits up 
to time t is simply the product Y\f^iP{xi,t). This is precisely the probability that the stopping time tg > t. Thus, 
the cumulative distribution of the stopping time tg, given the initial positions x, reads [31i. i39i] 

N 

(8) 



Prob[t.>t|f]=n-f(^) 



For large t, this cumulative distribution has a power law tail 

Pr„b[,.>,|f].(i)"'^£li|^^ (9) 

This asymptotic for survival probability tells us that the persistence exponent \% 9 — N /2. 

The result and its asymptotic counterpart © are exact. Next comes the heuristic part. We note that for 
large tg, the typical maximal displacement m in time tg must scale as: m ^ \/t^- Taking this relationship between 
the two random variables m and tg seriously, we see that Prob [m > L|a;] ~ Prob [tg > L'^\x] for large L. Since 
Qn{L\x) = Prob [m < , we conclude that for large L 

1 - Qn{L\x) = Prob[m > L\x] - Prob[t, > L^\x] - ^1^2^-'^^ (10) 
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where we used the result in Eq. ([9]) . Taking derivatives with respect to L and putting L 
behavior of the probability density P{m\x) of the maximum for large m 



PN{m\x) 



Xi X2 



. ,Xn 



m 



N+l 



m then gives an approximate 



(11) 



the 



as announced in (j4|). This heuristic scaling argument thus provides, up to an overall N dependent prefactor Bjq, 
leading asymptotic power law tail of the distribution of m in (|4|. 

To compute the prefactor i?jv exactly for any N one needs to go beyond this scaling argument. This requires a 
more sophisticated mathematical analysis that is carried out in the rest of the paper. However, it is possible to refine 
this heuristic argument, as shown below, that even provides the prefactor Bj^ exactly for large N . 

In the argument above, the main approximation was to replace m ^Jtl for large m and then use the exact 
asymptotic distribution of ts in Eq. Q to compute the tail of the distribution of m. This approximation clearly 
ignores fluctuations of m for a fixed tg. We now use an extreme value argument to show that this approximation 
actually becomes exact for large N . We consider again a group of Brownian motions starting at the initial positions 
X and examine their trajectories over a fixed time interval [0, tg], with m denoting their global maximum in [0, ts\. To 
compute the cumulative probability QN{L,ts\x) = Prob [m < L|af] , we consider the trajectories that stay below the 
level L till tg and also above the level 0. Now, for large L, the trajectories that contribute to Qn{L, ts\x) typically have 
large excursions. So, to a first approximation, one can ignore the lower boundary at 0. For the i-th walker, starting 
at Xi, the probability that its maximum stays below L can be computed easily: it is just the survival probability 
p{xo,ts) in Eq. ([T]) with the initial position xq = L — Xi. Thus the joint probability that all walkers stay below L till 
ts (ignoring the lower boundary at 0) is just the product 



N 



QNiL,ts\x) ~ J^crf 



L — Xi 



(12) 



For large argument, the error function behaves as 1 — erf(x) = e ^ /(x^/tt) ~ e ^ to leading order. Hence, for large 
L and large iV, one can write 



(L, ts\x) ~ exp [_iVe-^'/2*=l f [{L - aN)/b 



(13) 



where the scale factors = \/2ts In(A^) and = ^/T^JJ2AnJV) and the scaling function f{x) = exp[— e~^] is the 
standard Gumbel function. The pdf of m is just the derivative of the cumulative distribution. The derivative of the 
scaling function f'{x) = exp [—x — e^^] has a peak at a; = 0. This indicates that the random variable m has a peak 
at m = UN — \/2ts \n{N) and the width of m around its peak is = y/ts/i'^ IniV) that actually decreases with 
increasing N. Thus, for large iV, the random variable m approaches to the its mean value ajv with probability 1, i.e., 
m = y^2ts In(A^) with fluctuations around this value essentially negligible for large N. Using this relation in Eq. ^ 
provides the following tail for the cumulative distribution of m for large TV 



1 - QNiL\x) = Prob [m > X| f] ~ Prob 



1^2 

ts > — (ln7V)-^|f 



A 



Xi X2 



, Xn 



N- 



L 



N 



(14) 



with the prefactor An — [4 ln(iV)/7r]^^^. Taking derivative with respect to L then gives the tail of the pdf of m in 
Eq. with the prefactor, for large N, 



4 

-ln(iV) 

TT 



N/2 



En = NAn ^ N 

We will see later that the same asymptotic result also follows from a more rigorous approach. 



(15) 



III. A BACKWARD FOKKER-PLANCK METHOD: LAPLACE'S EQUATION 

In this section we show that the cumulative distribution Qn{L \ x) = Prob [m < L| a?] of the maximum m satisfles 
a Laplace's equation in an iV-dimensional hypercube with appropriate boundary conditions. To see this, it is flrst 
useful to consider Qfq{L\x) = Qn{x\L) as a function of the coordinates x for a given flxed L. Note that Qfq{x\L) is 
the probability that starting at x, the maximum of the process till ts stays below the level L. The idea is to derive 
a differential equation for Qfq{x\L) using a backward approach where one focuses on the evolution of the system via 



6 



([3]) over a small time interval [0,dt] starting from the initial positions x. According to ([3)), in this small time interval 
dt, the i-th particle moves from Xi to a:/ = Xi + r]i(0)dt where r]i(0) is the noise at t — that kicks the i-th particle. 
Now, starting from this 'new' initial coordinates the maximum of the system has to subsequently stay below L till 
the stopping time. Finally, one must sum over all possible values of the new coordinates Xi' . Thus one must have 

Qn [xi,X2, ■ ■ ■ :Xn\L) = {Qn {xi+ r]i{0)dt,X2 + r]2{0)dt, . . . ,xn + 7yAr(0)c?t|L)) (16) 

where the (.) denote the averages over the initial noises ?7j(0). Expanding the right hand side in a Taylor series and 
using (i) {r]i{0)) — and (ii) (?7i (0)77^(0)) — 1/dt (which follows from the delta correlator), one finds that Qn{x\L) 
satisfies the Laplace's equation in the A^-dimensional hypercube, < < L, 

\7^Qn{xuX2,...xn\L) ^0. (17) 

The information about the maximum is captured in the boundary conditions. For example, if Xj = for any i, 
Qat = 1 since if the i-th particle starts at the origin, the process stops immediately {tg = 0) and hence the maximum 
is necessarily (with probability 1) less than L. On the other hand, ii Xi ~ L for any i, one has Q = 0. This follows 
from the fact that if the i-th particle starts at L, it will immediately cross the level L and the probability that the 
maximum will stay below L till ts is necessarily zero. 

To summarize, Qn{x\L) satisfies Laplace'e equation ([T7]) in an A^-dimensional hypercube < Xi < L with Q — 1 
for any Xi = and Q = for any Xi = L. Thus, it reduces to an electrostatic problem where one needs to find the 
potential Qn{x\L) inside the hypercube [0,L]^, whose TV faces touching the origin are held at a constant potential 
Qat = 1 while the rest of the N faces are earthed {Qn = 0). We present the solutions explicitly for iV = 1 and N — 2 
in the next two subsections. 

Scaling: Let us remark that since the only length scales are the intial positions {xi} of the particles and the size of 
the box L, it is evident that the exit probability Qn{x\L) satisfies the scaling property 

Qn{x\L) = Qn (^^, Y'' ■ • ■ ' ^) ^ QNizi,z2, . . . , zn) (18) 
where the dimensionless scaled variables < Zi = Xi/ L < 1. 

Special initial condition and a duality relation: It is useful to consider a special initial condition where all the 
particles start from the same initial positions: Xi — x for all 1 < i < A^. In this case, the exit probabality is a function 
of a single scaled variable Q < z = x/ L <1: Qn{xi = x,X2 = x, . . . ,xn = x\L) = qN{x/L) where the scaling function 
qN(z) satisfies 

qN{z)+qN{l- z)^l. (19) 

This duality relation states that the exit probability through the right boundary at L starting from the initial position 
L ~ X (of all the particles) is exactly identical to the exit probability through the left boundary starting from the 
initial positions x (of all the particles). The duality relation (|19p in particular dictates that 

qN{z = 1/2) = 1/2 (20) 

for all N. The general solution of the Laplace's equation must satisfy Eq. (|19[) which actually provides a useful check 
for the validity of the solution. 



A. N = 1 



For N = 1, we have a second order ordinary differential equation 

d^Q 



dx\ 

whose solution is 



with Qi(0|L) = 1 and Qi(L L) = 0, (21) 



Qi{xi\L) = I - {)<xi<L (22) 
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Q2=i 



FIG. 3: Laplace's equation to be solved in a square (L x L) with boundary conditions as shown. 



Since Qi{xi\L) = Qi{L\xi) = Prob[TO < L|a;i], it follows, by differentiation that the probability density of the 
maximum m has a strict power law form for all m 



Pi(m\xi) 



for m > xi 



(23) 



which is normalized to unity over m £ [xi, oo] and all its integer moments diverge |35| . Thus the prefactor Bi = 1 in 

Clearly for = 1, the exit probability Qi{xi\L) is only a function of the scaled variable z — xi/L: Qi{xi\L) = 
qi{xi/L) where the scaling function qi{z) is simple 



qi{z) = 1 - z 



(24) 



and evidently it satisfies the duality relation (fT9|) . 



N 



The solution of the Laplace's equation becomes more involved in higher dimensions iV > 1. For N = 2, one needs 
to solve Laplace's equation in an (L x L) square with boundary conditions shown in Fig. [3l The solution can be 
written down explictly for this case, although the expression is rather cumbersome. Using separation of variables one 
finds the appropriate solution (46j 



i(xi,X2\L) = 1 



xi 2 \ ^ 



4)"-i [sinh(7i^(l-f)) + (-l)"sinh(^)] . fmrxi 



sinh(n7r) 



\ L J 



(25) 



It is straightforward to verify that ((25|) is a solution of the Laplace's equation. Next we need to check that it 
satisfies the 4 boundary conditions (see Fig. [S]). It is easy to check the two conditions: (i) Q2ixi = 0,a;2|i) = 1 for 
all X2 and (ii) Q2{xi = L,X2\L) = for all X2- The other two conditions can also be verified. For instance, putting 
X2 = L in (I^Sj) we get 



Q2{xi,L\L) = 1- >-sm — — 



(26) 



Using the identity 



1 . /mrxi\ TT / xi\ 



(27) 
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we verify that Q2{xi, L\L) = 0. Similarly, putting 2:2 = in (j25|) and using the identity 

n— 1 

one verifies the last boundary condition Q2{xi,X2 = 0|L) — 1. Note that even though the solution Q2{xi^X2\L) in 
(|25p is not manifestly symmetric under the exchange of Xi and X2, it is actually symmetric in xi and X2 as it should 
be. Later in Section 4, we will derive an alternative expression via the path integral method which is manifestly 
symmetric in xi and X2- 

Let us analyse the large L behavior of (j25p in the situation when xi and X2 are kept finite. Since Q2{xitX2\L) = 
Q2{xi/L, X2/L) is a function of only the scaled variables zi — xi/L and Z2 — X2/L, the i — >■ 00 limit is equivalent to 
taking the limits zi and Z2 0. Clearly Q2{zi = 0, 2:2 = 0) = 1. It is easy to check that the first derivatives dz^Q 
and dz2Q vanish at the origin (zi = 0, Z2 = 0). Similarly, the second derivatives dl_^Q2 and dl^Q2 also vanish at the 
origin. So, in a Taylor expansion of (32(2^1,^2) the first nonzero term is the cross derivative, indicating the following 
leading order behavior as zi — >■ and Z2 — )■ 

Q2{ZI,Z2) =1- A2Z1Z2 + ... (29) 

The amplitude A2 is given by 

- 1 = /TT ^ \ — ly n corn^r^7^J — , 

^-^ ^-^ sinhfnTr) 

.n=l n=l ^ 



dzidz2 U=o.^2=o 



(30) 



where the last line follows from (j25p . Fortunately, the sums can be carried out explicitly using some known identi- 
ties [131 to yield 



A2=g^[r(l/4)]^ = 2.18844... (31) 



Hence for large L 



Prob [m < L|a;i, 2:2] = Q2{xi,X2\L) ~ 1 - A2 ^^^^ . (32) 

which leads to the announced power law tail for the probability density of the maximum m 

P2{m\xx,X2) B2^^^^\ with B2 = 2A2 = 4.37688 .. . (33) 

To compare with the TV = 1 case, let us consider the special initial condition where both particles start from the 
same initial position: x\ — X2 — x. In this case, the exit probability Q2{x\ — x,X2 = x\L) = q2{x/L) with scaling 
function 92 (^) given by 

, 2 ^ (-1)"-Msinh(n7r(l - z)) + (-l)"sinh(n7rz)] . , 

q2{z)^l- z+ -y ^ ^- ^ ^ / / ^ sin nTTZ . 34 

TT ^-^ n smh TiTT 

n— 1 ^ ^ 

One can verify that (72 (^) satisfies the duality relation ([T9l) . Near z — and z = 1, the scaling function has the 
asymptotics 



I 1 - A2z2 as z ^- 0, 
'''^'^^^A2{l-zf asz^L 



Comparing with the expression for the = 1 case, Eq. (p4)) . one finds that q2{z) has a much richer functional 
form. An analytic prediction for the function q2{z) is in excellent agreement with results obtained from Monte Carlo 
simulations (Fig. [5]). 



IV. EXACT SOLUTION FOR ALL BY PATH COUNTING METHOD 



The approach based on the Laplace equation (Sect. IIII[) is difficult to extend for > 2. The technical problem is 
to find the exact solution of the Laplace equation V'^(3Ar(zi, Z2, . . . , z^r) = in the TV-dimensional hypercube of the 
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scaled variables < = x^/L < 1, with the prescribed boundary conditions that Qn = 1 for any Zi = and Qjv = 
for any zi = 1. For A'^ > 2, it is not easy to find an explicit solution to this problem. In this section, we use an 
alternative path counting method that is physically more explicit and, in addition, allows us to write down the exact 
solution QAr(zi, Z2, . . . , zn) for all N. This method thus provides an alternative way to solving the Laplace's equation 
in a hypercube with the prescribed boundary conditions. 

To set up the path counting method for general N , we need two basic ingredients from the single particle problem 
in a box [0, L] with absorbing boundary conditions at the two boundaries and L: 

(i) the survival probability S{xo,t,L) which counts the probability that a single particle, starting at xq at time 0, 
remains inside the box [0,L] during the time interval (0,t). 

(ii) the first passage probability density F{xo,t, L) that denotes the probability density that the particle, starting 
initially at xq {0 < xq < L), exits for the first time the box through the boundary at (and not through the other 
boundary at L) at time t. 

One can write the survival probability S{xo,t,L) — G{x,XQ,t, L) dx where G{x,xo,t, L) denotes the Green's 
function counting the probability density that the particle reaches x at time t, starting from xq at time 0, while 
staying inside the box [0, L] during time t. This Green's function can be computed (a) either by solving the diffusion 
equation dtG — Dd^G (the diffusion constant is = 1/2 for our choice of the noise term) with absorbing boundary 
condition G{x = 0,xo, t, L) = G{x = L, xq, t, L) = and the initial condition G{x, xq, 0, L) = S{x — xq), or (b) by the 
path integral method. The resulting Green's function admits the following representation 

Gix,xoAL)^l±sin{;^) sm(!^) e"^*. (36) 

ji=i 

The survival probability, after integrating over the final position x, is 

Q( , T\ 2 ^ [1 - (~1)"] . fn7TXo\ 

S{xo,t,L)^-} sm — — e . 37) 

n=l 

To calculate the first passage probability density F{xQ,t) through at time t, we note that this just counts the 
fiux of particles going out of the box through at time t. The flux through a point x in the positive direction is the 
probability current —DdxG (with D = 1/2). Hence, the first-passage probability through the origin that counts the 
fiux through the origin in the negative direction is simply, F{x(),t) — DdxG\^_^. Using G from we get 

CO 2 2 

F{xo,t) = ^^Y.^sin['^) e-^\ (38) 

n=l 

Armed with these two ingredients from the single particle problem, we are now ready to compute the exit probability 
Qn{x\L) for the A^-particle problem. Consider first the event shown in Fig. [5] where one of the particles (say the 
i-th one) exits the box for the first time between time tg and tg + dtg while the — 1 other particles stay inside 
[0,i] till this time tg. Clearly, the probability for this event, using independence of walkers, is given by the product 
F{xi, tg) Yij^i S{xj,ts) dtg. Now, the particle that hits (whose label is i) can, in fact, be any one of the N particles. 
Hence we have to sum over the hitting index i from 1 to N. In addition, this event can occur at any time tg, so we 
need to integrate over tg . This path counting method then gives the following net contribution to the exit probability 

Qn{x\L)^J2 F{x^,tg)Y[S{x,,tg)dtg. (39) 

Substituting the results for S{x,t) and F{x,t) respectively from Eqs. ([57| and integrating over tg and setting 
Zi = Xi/L gives our main exact result, valid for all N, 

r, f^\T\ r, f ^ f ST^ a(ni, ^2, . . . , tiat) sm{nkTrzk) , 

Qn{x\L) = Qn{zi,Z2,...,zn) = [- ] V — ; — r (40 

^'"^ + + ^^^^^ nk 

where each index = 1, 2, . . . for all 1 < i < iV and 

N 

a{n,,n2, . . . , oat) = ^ [1 - (-1)"^] . (41) 
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Equation (l40l) is an exact solution of the Laplace's equation in the iV-dimensional hybercube < < 1 with the 
prescribed boundary conditions which has a virtue of being manifestly symmetric with respect to the interchange of 
the coordinates Zi's. 

To extract the behavior in the limit of large L, with x^'s fixed, we need to take the limit — > in Eq. PO)) . This 
limiting behavior reads 

Qn{zi,Z2, . . . , zn) = 1 - AnZi Z2...ZN + 0{zf) (42) 



The coefficient An can be computed by taking the derivatives of Eq. (|40| with respect to each Zi and then putting 
Zi = for all i. This gives 

An = -A^2^ lim „2 ,^2: 1'2 n [1 - (-l)"1cos(n,.z,). (43) 

Noting that (3Ar(af|i) = Qn{zi, Z2, . . . , zat) is precisely the cumulative distribution of the maximum m, one thus 
gets 

Prob[m <L\x]~l- An ^^^^j^^'^^ asL^oo. (44) 

Differentiating with respect to L and putting L ~ m gives the exact power-law tail of the pdf of the maximum m in 
Eq. (j4|) with the prefactor Bn is given by the formal sum 



Bn = NAn = -N' 2^ lim Y: "i jl"^""^) Jl [1 _ (_l)".] cos(n,.z,) (45) 



where each index rii runs over all positive integers. 

The formal sum in Eq. (|45p can be explicitly evaluated for = 1 and N = 2. For = 1 



Bi = -2 lim V cos(ni7rzi) = 1 (46) 

21^0 ^ 

711 = 1,2,... 



in agreement with already known results. For N = 2, equation (|45|) gives 

„ , ■^-^ cos(ni7rzi) , / s„ , , , 1 

B2 = -16 lim y ^ 2 2 1 - -11 cos n27rz2 = 



n 4 



= 4.37688 . . . (47) 



In Appendix A, we show how to compute the above sum explicitly. 

For N > 2, we haven't deduced explicit expressions for the sum in (145^ . However, one can reduce it to a simpler 
form where the sum is rapidly convergent and can then be evaluated by Mathematica. As an example, for = 3, 

9 4, , ^ n2 + n2(-l)"2 

B3 = -6 + -^r'*(l/4) + 727r 2J ^^)-^-==:^ = 15.3369 . . . (48) 

^'^ m,n2 V"? + nl sinh f Trv^nf + 7 



In the limit of large iV, one can evaluate the formal sum (see Appendix B) to obtain the limiting behavior 

Bm^N 



- In(iV) 

TT 



N/2 

(49) 



in perfect agreement with the heuristic result in Eq. ([15 



Special initial condition: Finally, let us consider the special initial condition when all the particles start from the 
same point: xi = X2 = ■ ■ ■ ~ x where < a; < L. In this case, the exit probability in Eq. (j40p reduces to a function 
of one scaled variable z = x/L: Qn{x, x, . . . ,x\L) = q^iz) given by 



/^ AT f '^\^ Sr^ nisin(ni7rz) -p-r sin(nj7rz) 
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One can check that qN{z) satisfies the duaUty relation qN{z) + qni^ — z) = \ and thus (7Ar(l/2) = 1/2 for all N . 

Performing this multiple sum directly by Mathematica is difficult as it converges slowly. To circumvent this problem, 
we first perform the sum over ni in Eq. ([50| using the following identity ^\ 



E 

k=l 



k sm.{k'Kx) 



TT sinh {^{x — a)) 
2 sinh(7ra) 



(51) 



This gives 



TT 



N~l 



E 



sinh 



7r(l - z).Jnl + nl + 



712 ,71^,. ..,n]^-^odd 



sinh TTy'' 



"N 



N 

n 



sin(rtj7rz) 



(52) 



The multiple sum in Eq. (|52l) is now rapidly convergent and can be easily evaluated by Mathematica. In Fig. [2] we 
plot this function qN{z) for TV = 2 and iV = 3. For = 2, it of course coincides with the earlier expression ([M)) 
obtained via the Laplace's method in the previous section. 

Average Maximum: For fixed identical initial position of all particles Xi = x, another interesting question is how 
does the average maximum (till the stopping time ts) depend on A''? For iV = 1 the average is infinite but for all 
> 1 it is finite. However, does the average maximum for A'^ > 1 increase or decrease as the number of walkers 
A" increases? The answer to this question is not intuitively obvious. However, knowing the function q^iz), one can 
compute the average maximum in the following way. We have, Prob[m < L\x] = qN{x/L) where qN{z) is given in 
Eq. ([5^ . Therefore the probability density function of m reads 



PN{m\x) 



^'^'^ (^) ' 



m > X 



where q'j\j{z) = dqi^{z) / dz. The first moment is then given by, for all x > and A^ > 1 



/•oo ^ / X \ 

I — q'pf ( — ) dm ^ Cn X 



with prefactor 



Cn — — 



dz. 



(53) 



(54) 



(55) 



Thus the average maximum, for A^ > 1, is proportional to x for all x and the proportionality constant Cjv is given by 
Eq. (ISSl) . We were unable to carry out the integral in Eq. (ISSl) in closed form. However, it is clear that as A^ — >■ oo, 
using q^iz) — 0{l/2 — z), one gets Cn — >■ 2. On the other hand Cjv diverges as A^ — ^ 1. Thus, Cjv decreases when A^ 
increases. These results are supported by Monte Carlo simulations (Fig. [5]). 



V. NUMERICAL SIMULATIONS 



Our analytical result for the exit probability through the origin q]\[[z), when all particles start at the same scaled 
position < z = x/L < 1, is tested by Monte Carlo simulations. For a fixed box of size L, the exit probility for the 
A'^-particle problem can be efficiently computed using a method proposed in Ref. [i^. Naively, to compute q]\f(z), one 
would first fix the starting point < z < 1 of all the Brownian motions and then generate different realizations of 
the process and compute the fraction of realizations where the first exit happens through the origin. One would then 
repeat the procedure for each value of < z < 1 (with an appropriate bin size) in order to compute the full function 
qN{z) over the range z e [0, 1]. Instead of repeating the simulation for each starting point z, it turns out to be more 
efficient to follow a different algorithm described briefly below. 

For illustration, we take the example of just one Brownian motion (A^ = 1). The method is easily generalized for 
all N. We start the Brownian motion at the origin, let it evolve in time and record the maximal (xmax) and the 
minimal (a^min) position reached by the walker up to time t. The process is halted when a;inax — a^min > L for the first 
time (see Fig. 14]). Keeping a^min and Xmax fixed, we now horizontally slide the whole configuration thereby changing 
the starting point (see the second and the third panel of Fig. 01). Measuring all distances with respect to a;niin, it is 
then clear that this configuration contributes 1 to qi{z) for z € [0, 1 — Xmax/i] and for z G [1 — Xmax/ L,l]. So, 
for this configuration, we just record the number y = 1 — Xmax/i and the fact that it contributes 0[y — z] to qi{z). 
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We then repeat this procedure for another configuration starting at 0, do the shding and record the value of y and 
the associated 9{y — z). We repeat this, say, for Ng number of samples. To sum all the contributions, we first sort 
the values of y associated with the configurations in increasing order. Let {yi,y2, ■ • ■ denote the sorted values. 

Then, we assign (?i(yi) = 1, qi{y2) = 1 ^ ^/Ng and in general qi{yi) = 1 — {i — 1)/Ns. This generates the full curve 

Thus this algorithm has two advantages: (i) it does not require to repeat the simulations for each value of z and 
(ii) it does not require any specific choice of bin sizes. This clearly makes the simulation much faster. The results of 
our simulations for = 2, = 3 and = 6 are plotted in Fig. ([5]). The agreement with the analytical prediction is 
excellent. For all N, qN{z) satisfies the duality relation: qN{z) + 9^(1 — z) ~ 1. 

In Ref. [i^ we had studied the exit probability from a box [0, L] through the upper boundary L (rather than the 
lower boundary considered here) for a single particle whose motion is described by a generic self-affine stochastic 
process. This probability was called the 'hitting probability' in [isj . In the present paper, we are considering the 
complimentary event of exiting the box [0,L] through the lower boundary 0. In the notation of the present paper, 
the hitting probability (of the boundary L) would correspond to 1 — qN{z) HOl- For a generic self-affine process, it 
was shown in Ref. |43] that the exit probability should have a power law behavior close to the origin 1 — q{z) ~ z"^ 
with (f> ~ 9/H as mentioned in the introduction. In addition, it was observed in ref. |43| that for many processes 
(but not all), once we know the exponent </>, the full function q{z) over the range z e [0, 1] is described by a universal 
one-parameter (parametrized by (f)) form (43| 

q^z) = 1 - 4(0, 0) = 1 - M ^'[^,(1 _ u)]-l'-'du. (56) 

The function 1^(4', (f>) is the incomplete regularized beta function. In our present problem, we have seen in Eq. (|42p 
that for small z, qN{z) 1 — Anz'^ indicating <j) = N . It is then natural to investigate if our result for the full 
function qN{z) can be re-expressed as the universal functional form in Eq. (|56| with cf) — N. Interestingly, the answer 
is no, as it is clearly shown in Fig[6l thus providing us with a counterexample. 

We have also computed the prefactor numerically from the Monte Carlo simulations up to = 7. The results 
are shown in Fig. [7] by squares. For A^ = 1, A^ = 2 and A^ = 3, they agree with our exact analytical predictions 
Bi = 1, i?2 = 4.37688 . . . and B^ — 15.3369 ... as discussed in Section IV. Had the qN{z) be described by the universal 
function in Eq. ([56|) with (f) — N, one would get from the small z expansion in Eq. (fSH)) . qN{z) -> 1 — A^z^ with 
An = T{2N)/NT'^{N). This would predict the prefactor Bn = NAn = r{2N)/r'^{N). In Fig. H this prediction from 
the universal curve is shown by the dashed line. Clearly, it does not match the simulation results, confirming once 
more that qN{z) is not described by the universal function in Eq. (I56p . We have also plotted the exact asymptotic 
prediction oi Bn in Eq. (|49| as a solid line for comparison. While it is difficult to extract the small z behavior of 
qN^z) and hence Bn for larger values of A^, we notice that the asymptotic large A^ behavior is already approached 
for A^ = 7. 

Finally we have also computed the average maximum (till the stopping time ig) of Af Brownian motions starting 
from the same initial positions a; > 0. We verified that for all A^ > 1, the average maximum exists, and it is 
proportional to x, as predicted analytically in Eq. (|54p . In Fig. [5] we plot the amplitude Cjv, given analytically in 
Eq. (|55|) . as a function of A^. As predicted, Cn approaches 2 for large A^. 
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FIG. 5: The exit probability qN{z) through the origin obtained from Monte Carlo simulations for N = 2 (red squares), = 3 
(deep blue circles) and N — 8 (blue triangles). The numerical data are compared to analytical results (solid lines) for N = 2 
and A'' = 3 (obtained from Eq. (|52|) respectively for N = 2 and = 3). 
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FIG. 6: Difference between the exit probability for N = 2 and N = 3 and the universal function in Eq. (I56|l with <f) = N. Solid 
lines are the analytical results from Eq. (|52|) with N — 2 and TV = 3. Symbols are the numerical results for 10"" realizations of 
N = 2 (squares) and A*' = 3 (circles) Brownian motions starting from the same initial position in a box of size L = 1000. 



VI. CONCLUSION 



To summarize, we have presented an exact solution for the probability distribution of the maximum m of a set 
of N independent Brownian motions starting at the initial positions x = {xi,X2, ■ ■ ■ ,xpf} on the positive half-axis 
and the process terminating when any one of the walkers crosses the origin. We have shown that for large m, the 
pdf of m decays as a power law, Ppf{m\x) ~ Bn{xiX2 ■ ■ . xat)/™^"*"^ where the prefactor Bn has an interesting N 
dependence. For a fixed > 1, integer moments of m up to order {N — 1) are finite, while all higher integer moments 
are infinite. The cumulative distribution of this maximum also provides an exact solution to the first-exit probability 
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FIG. 7: Dashed line: Bn = r(2A'^)/r (A'^) from the universal function in Eq. (|56p . Solid line: aymptotic behavior for large 
N given in Ea. (l49|) and squares: direct simulations of A'' Brownian motions starting from the same position in a box of size 
L = 1000. Averages are performed over l(f samples. 




FIG. 8: Monte Carlo results for the average maximum of Brownian motions till the first passage time through the origin, all 
starting at the same initial position x > 0. The coefficient Cjv ~ {m)/x, see Eq. (|55p . is plotted vs. A'^. Averages are performed 
over 10^ realizations. 



through the origin (rather than through L) of N walkers from a box [0, -Zv]. Incidentally, our path counting method 
also provides an exact solution to the iV-dimensional Laplace's equation V'^Qn = in a hypercube [0,L]^ with the 
boundary conditions Qat = 1 on any face of the hypercube passing through the origin and Qat = on the rest of the 
faces. Monte Carlo simulations confirm our analytical results. 

This work raises some interesting open questions. We have focused only on the maximum m till the stopping time 
ts of TV independent walkers. Another interesting observable is not just the actual value of the maximum, but the time 
tm at which this maximum occurs before the stopping time t^. This random variable has recently been studied in a 
number of contexts. For a stochastic process over a fixed time interval [0,t], the distribution of the time tm has been 
computed for a variety of Brownian paths, such as a free Brownian motion, Brownian bridges, Brownian excursions 
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and Brownian meanders, using path integral method ^ and also by an alternative functional renormalization group 
method J^j. The distribution of tm was also computed exactly for the random acceleration process which is a non- 
Markov process fsil. It has also been computed both for independent Brownian walkers [H, and very recently 
for vicious walkers '5z]. On the other hand, when the process stops at a random stopping time tg where for instance 
ts is the first time a walker hits the origin, the distribution of tm has been computed for a single Brownian motion 
{N — 1) where it is already nontrivial [53|. It would be interesting to extend the results of Ref. [53] to the case of 
> 1 independent Brownian motions. 

Another interesting challenging problem would be to compute the distribution of m as well as that of i,„ for a set 
of vicious walkers till the stopping time ts when the walker closest to the origin crosses the origin for the first time. 
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Appendix A: Evaluation of B2 

To evaluate the sum in Eq. ()47|) we first split the sum into two parts: B2 = Ii + 12 with 

Ii = —16 lim > — i — 7j cosfriiTTZi) cos(n27rz2) (Al) 
zi,z2->0 ^-^ ni + 

— ^ ^ cosfriiTTZi) cos(rt27rz2) (A2) 

^1.^2 ni + ni 

ni,n2 

where all the sums run over positive integers. 

Let us first evaluate Ii . Due to the symmetry we can re- write Ii as 



Ii ~ —8 lim \^ — i 1 cosfniTTZi) cos(n27rz2) ~ ~8 lim } cosfriiTrzi) cos(n27rz2) = ^2 

2i,z2^o ^-^ ni + ni zi,z2^o 



(A3) 



In the last step we have used the identity, X]n>i cos(n7rz) — —1/2. 
To evaluate I2 we use another standard identity [l^l 

E°° (— 1)*^ cos(fc7rz) TT cosh(a7rz) 1 
k^ + a^ " 2^ sinh(a7r) " 2^" ^ 

k—l 

We now sum over n2 in Eq. (jA2p using above identity. This gives 

CX3 00 00 

/2 = — 8 lim cos(ni7r2i) -f 87r hm — cos(ni7rzi) cos(ni7rz2) = 4 + 87r — -. (A5) 

zi^o ^ zi,z2^o ^ smh(ni7r) ^ smh(ni7r) 

ni=l ni = l ^ ' ni = l ^ ' 

The remaining sum in (|A5P can be explicitly evaluated using the identity [l^l 

V " =J_T^(l\^l. (A6) 
^ sinh(n7r) 327r3 1 4 / 47r ' ^ ' 

Adding Ii and I2 we arrive at the announced expression ([5]) for 82- 



Appendix B: Evaluation of Bm for large TV 



Here we show that in Eq. (|45|) behaves, to leading order for large N, as 



B 



N 



N 



- In(iV) 

TT 



N/2 



(Bl) 
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First we rewrite the denominator in Eq. (j45p using the integral representation 

1 



n? + n2 + . . . + n% 



Using this representation, one can then decouple the sums over different indices in Eq. (1451) giving 
Bn = -A^^ 2^ lim / dt "^nj cos(ni7rzi) e~"i* H XI ~ (-1)"'] cos(nj7rZj) e"*": 

-^0 ni=l 7#ln,=l 



(B2) 



(B3) 



If we now take the limits 



inside the sums, each of the sums is convergent. However, it is easy to check that the 



integrand, as a function of t, diverges as ~ /;i+JV/2 f ^ q Thus the integral is longer convergent. To circumvent 
this difhculty, one can use a standard regularization scheme used often in evaluating the Madelung constant in the 
context of lattice sums in crystals [4§| whereby we introduce a parameter s and rewrite Eq. (jB3P as, upon taking the 



limits 



inside the sums. 



Ar2 oiV foo ^ 



N-l 



(B4) 



Note that this integral is covergent for all s > (1 + N/2). The idea is to first evaluate Bn{s) for large N with 
s > (1 + N/2) and then analytically continue this result to s — >■ 1 to evaluate Bn — Bn{s — >■ 1). 
Let us next define the function 



,-(2fe+l)^t 



(B5) 



k=0 



Up to the factor \/t/n, this function g{t) can be expressed in terms of standard Jacobi theta functions [4^. Then 
one can rewrite Eq. (jB4p as 



2N^{V^) 



N-l /-oo 



r(.) 



E 



(B6) 



.71=1 

To evaluate this integral for large N, we need to know how the function [g{t)]^^^ behaves for large N. 

Let us first focus on the function g{t) in Eq. (jBSp . Clearly, for large t, the dominant contribution comes from the 
k = term in the sum and hence g{t) ^ Vie^* as t — >■ c». In contrast, the opposite limit f — ^ is more tricky. To 
derive its behavior as i — >■ 0, we first use the following Jacobi identity [49l | 



oo j — oo 

1 + 2 ^(-l)'=e-'=' - = 2 e-(''=+i)'"' 

fe=l ' fc=0 

to rewrite the function g{t) (upon identifying t — 7r^/4z) as 

oo 

g{t) = l + 2^(-l)'=e-"''='/4t^ 



/4z 



(B7) 



(B8) 



fe=i 



This representation of g{t) is more amenable to the asymptotic analysis in the i — > limit. We obtain from Eq. 
the leading behavior of g{t) as i -> 
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(t)~l-2e-'^'/4* + 



(B9) 



Thus g{t) has an essential singularity at t = and it approaches to 1 as t —> in an extremely flat way. Thus the 
function g{t) starts at g{0) = 1, stays flat for a while and then decreases exponentially fast to as i increases (see 
Fig. ®). 

Next consider the function [g(i)]^ that appears in the integral in Eq. (|B6p . As N increases, the function [^(t)]^ 
almost approaches a step function (see Fig. (©) 



[g(t)r ^0{t*{N)-t) 



(BIO) 
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FIG. 9: The function plotted as a function of f for iV = 1 (black solid), N = 5 (red dotted) and iV = 50 (blue dashed). 

As N becomes large, [g{t)]^ almost approaches to a step function. 

where the characteristic scale t*{N) decreases very slowly with increasing N. One can easily estimate t*(N) for large 
N from the asymptotic behavior in Eq. (jB9l) . For small t, one finds, to leading order for large N 



N 



l-2e-^'/4* 



N 



exp 



(Bll) 



Thus as increases, it approaches to rapidly for all t > t* where 2Ne sa 1. This provides an estmaite of 

t*{N) which, to leading order for large N, reads 



41n(7V)' 



(B12) 



Therefore, for large N, using Eq. (|B10[) . we can cut off the upper limit of the integral in Eq. (|B6p at t = t*{N) 
and replace [g{t)]^~^ by 1 over the interval t S [0,t*{N)]. Furthermore, over this small interval t G [0,t*{N)], one 
can replace the function X^riLi e " * by its small t behavior ~ -JH j ^^i?!"^ . Substituting these results in Eq. (|B6P 
then yields, for large iV, 



B 



N 



^' 2r(s) Jo 2r(s) (s-iV/2-1)- 



(B13) 



In deriving this result, we have assumed s > (1 + N/2). After obtaining this large N formula for i?Ar(s), we can now 
analytically continue it to s — > 1 which finally yields 

Bn = Bn{s ^ 1) = [t*(7v)]-^/2. (B14) 

Upon using the expression for t*{N) from Eq. (jB12|) gives the final large N expression for Bpf in Eq. ()Bip . 



[1] K.L. Chung, Bull. Amer. Math. Soc. 81, 742 (1975). 

[2] D.P. Kennnedy, J. Appl. Probab. 13, 371 (1976). 

[3] R.T. Durrett and D.L. Iglehart, Annals of Probab. 5, 130 (1977). 

[4] D.A. Darling, Annals of Probab. 11, 803 (1983). 

[5] J. Pitman and M. Yor, Elec. J. Probab. 4, 15 (1999). 



18 



[6] For a brief review see S.R. Finch, Variants of Brownian Motion available at http:/ /algo.inria.fr/csolve/br.pdf 
[7] S.N. Majumdar, J. Randon- Furling, M.J. Kearney, and M. Yor, J. Phys. A. Math. Theor. 41, 365005 (2008). 
[8] A.N. Kolmogorov, Annals of Math. Statist. 12, 461 (1941). 
[9] N.V. Smirnov, Mat. Sbornik 6, 3 (1939) (In Russian). 
[10] F. PoUaczek, Comptes Rendus 234: 2334 (1952). 

[11] F. Spitzer, Trans. Am. Math. Soc. 82: 323 (1956); Duke Math. J. 24: 327-343 (1957). 

[12] E.G. Coffman, P. Flajolet, L. Flato, and M. Hofri, Probability in Engineering and Informational Sciences, 12: 373-386 
(1998). 

[13] A. Comtet and S.N. Majumdar, J. Stat. Mech.: Th. and Exp. P06013, (2005); S.N. Majumdar, A. Comtet, and R.M. 
Ziff, J. Stat. Phys. 122, 833 (2006). For a review see, S.N. Majumdar, "Universal First-passage Properties of Discrete-time 
Random Walks and Levy Flights on a Line: Statistics of the Global Maximum and Records" (Lecture notes for the Leuven 
summer school FPSP-XII (2009), arXiv: 0912:2586, to appear in Physica-A (2010), doi: 10. 1016/j.physa.2010. 01.021). 

[14] R.M. Ziff, S.N. Majumdar and A. Comtet, J. Phys. C: Cond. Matter 19, 065102 (2007); R.M. Ziff, S.N. Majumdar and A. 
Comtet, J. Chem. Phys. 130, 204104 (2009). 

[15] S. Raychowdhuri, M. Cranston, C. Przybyla, and Y. Shapir, Phys. Rev. Lett. 87, 136101 (2001). 

[16] S.N. Majumdar and A. Comtet, Phys. Rev. Lett. 92, 225501 (2004); J. Stat. Phys. 119, 777 (2005). 

[17] G. Schehr and S.N. Majumdar, Phys. Rev. E 73, 056103 (2006); G. Gyorgyi et. al. Phys. Rev. E 75, 021123 (2007); T.W. 
Burkhardt et. al. Phys. Rev. E 76, 041119 (2007); J. Rambeau and G. Schehr, J. Stat. Mech.: Th. and Exp. P09994 
(2009). 

[18] G. Schehr and P. Le Doussal, J. Stat. Mech.: Th. and Exp. P01009 (2010). 
[19] S. Carmi, L. Turg eman, and E. Bar kai, arXiv: 1004.0943 , 
[20] V. Tejedor et. al, larXiv:1001.4412l 

[21] G. Schehr, S.N. Majumdar, A. Comtet, and J. Randon- Furling, Phys. Rev. Lett. 101, 150601 (2008). 

[22] M. Katori, M. Izumi, and N. Kobayashi, J. Stat. Phys. 131, 1067 (2008); N. Kobayashi, M. Izumi, and M. Katori, Phys. 

Rev. E 78, 051102 (2008). 
[23] T. Feierl, Proc. of IOWA 2009, Lecture Notes m Computer Science, vol. 5874 (2009). 
[24] C. Nadal and S.N. Majumdar, Phys. Rev. E 79, 061117 (2009). 

[25] K. Johansson, Probab. Theory Relat. Fields 123, 225 (2002); M. Katori and H. Tanemura, J. Math. Phys. 45, 3058 (2004); 

P. Ferrari, Lecture Notes of Beg-Rohu Summer School, available at http://ipht.cea.fr/Meetings/BehRohu2008/ 
[26] J. Randon-Furling, S.N. Majumdar and A. Comtet, Phys. Rev. Lett. 103, 140602 (2009). 
[27] S.N. Majumdar, A. Comtet and J. Randon-Furling, J. Stat. Phys. 138, 955 (2010). 

[28] for a brief review on "extreme value statistics of correlated random variables", see S.N. Majumdar and P.L. Krapivsky, 
Physica A, 318, 161 (2003). 

[29] S.N. Majumdar, Brownian Functionals tn Physics and Computer Science, Curr. Sci. 89, 2076 (2005) (also available at 

arXiv: /cond-mat/0510064). 
[30] M. Tachiya, Radiat. Phys. Chem. 21, 167 (1983). 

[31] A. Blumen, G. Zumofen, and J. Klafter, Phys. Rev. B. 30, 5379 (1984). 

[32] S.F. Burlatsky and A. A. Ovchinikov, Sov. Phys. JETP 65, 908 (1987). 

[33] S. Redner and P.L. Krapivsky, Am. J. Phys. 67, 1277 (1999). 

[34] D. Dhar and R. Ramaswamy, Phys. Rev. Lett. 63, 1659 (1989). 

[35] M.J. Kearney and S.N. Majumdar, J. Phys. A. Math. Gen. 38, 4097 (2005). 

[36] M.J. Kearney, J. Phys. A. Math. Gen. 37, 8421 (2004). 

[37] W. Feller, An Introduction to Probability Theory and its Applications (New York, Wiley, 1968). 
[38] M. Bramson and J.L. Lebowitz, Phys. Rev. Lett. 61, 2397 (1988). 

[39] A.J. Bray and R.A. Blythe, Phys. Rev. Lett. 89, 150601 (2002); R.A. Blythe and A.J. Bray, Phys. Rev. E 67, 041101 
(2003). 

[40] G. Oshanin, O. Benichou, M. Coppey, and M. Moreau, Phys. Rev. E 66, 060101 (2002). 

[41] A.J. Bray, S.N. Majumdar, and R.A. Blythe, Phys. Rev. E 67, 060102 (2003); M. Moreau et. al., Phys. Rev. E 67, 045104 
(2003). 

[42] S.B. Yuste et. al. Phys. Rev. E 78, 021105 (2008). 

[43] S.N. Majumdar, A. Rosso, and A. Zoia, Phys. Rev. Lett. 104, 020602 (2010). 

[44] S.N. Majumdar, Persistence in Nonequilibrium Systems, Curr. Sci. 77, 370 (1999) (also available at 

|arXiv:cond-mat/990740 7 ) . 
[45] S. Redner, A Guide to First-Passage Processes (Cambridge University Press, Cambridge, 2001). 

[46] N.N. Lebedev, LP. Skalskaya, and Y.S. Uflyand, Problems of Mathematical Physics (Prentice-Hall, Englewood Cliffs, N.J, 
1965), 

[47] A. P. Prudnikov, Yu.A. Brychkov, and O.I. Marichev, Integrals and Series, Vol. 1 (Overseas Publishers Association, Ams- 
terdam, 1986). 

[48] D. Borwein, J.M. Borwein, and K.F. Taylor, J. Math. Phys. 26, 2999 (1985). 

[49] I.S. Gradshteyn, I.M. Ryzhik, Tables of Integrals, Series, and Products (Academic, New York, (1980)). 

[50] Please note that in ref. [43| the notation q{z) was used for the hitting probability (for exiting through L), but here we use 

qN{z) to denote the complimentary event of exiting the box through 0. We apologise if this causes any confusion. 
[51] S.N. Majumdar, A. Rosso, and A. Zoia, J. Phys. A: Math. Theor. 43, 115001 (2010). 
[52] J. Rambeau and G. Schehr, arXiv:1004:0141. 



19 

[53] J. Randon-Furling and S.N. Majumdar, J. Stat. Mech.: Th. and Exp. P10008 (2007). 



